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LIST OF SYMBOLS 



B m Viscous friction coefficient of the motor (oz-in/rad per sec) 

B]_ Viscous friction coefficient of the load (oz-in/rad per sec) 

Eba Back EMF of phase A (volts) 

E B b Back EMF of phase B (volts) 

Ebc Back EMF of phase C (volts) 

I A Phase A current (amperes) 

I 3 Phase B current (amperes) 

Iq Phase C current (amperes) 

1^ Power supply current (amperes) 

J M Moment of inertia of motor shaft (oz-in-sec 2 ) 

J;L Moment of inertia of load shaft (oz-in-sec 2 ) 
kb Back EMF constant (volts/rad per sec) 
k^, Torque constant (oz-in/ampere) 

L Inductance of the stator winding (Henrys) 

Iab* ^CA’ ^BC Loop currents (amperes) 
k a dj Air-gap flux adjustment factor 

k„i Magnitude of speed versus torque curve slope (rad per sec/oz-in) 
R a Resistance of the stator windings (ohms) 

R s Power supply interval resistance (ohms) 

4>A » 4>b» 4>c Per P hase air-gap flux (webers) 

Rq 'On' resistance of the power transistor (Ohms) 

T m Motor restraining torque (oz-in) 

T Load restraining torque (oz-in) 

T A . Tb, Tq Per phase developed torque (oz-in) 
w m Speed of the motor shaft (rad/sec) 
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GR Gear ratio 



6 Angular displacement (radians) 

V A , Vg, V C Phase terminal voltage (volts) 

V5 Power supply terminal voltage (volts) 

RPSA , RPSB, RPSC Position sensors logic signals 
PWM Pulse width modulation 
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INTRODUCTION 



Recent improvements in rare-earth magnetic materials for use in brushless 
dc motors have allowed reconsideration of electro-mechanical actuator systems 
for applications requiring very high ratios of torque-to-inertia . The 
investigation discussed herein has been concerned with characterizing mathe- 
matically the dynamical features of a missile fin actuation system, from the 
input to the brushless dc motor to the output shaft of the mechanical actuator. 
The physical model is based upon an existing prototype actuator currently 
under evaluation at the Naval Weapons Center, China Lake, California. 

In general, brushless dc motors produce torque through the interaction of 
a magnetic field generated by a permanent magnet rotor and a dc generated 
magnetic field in the stator. The rotating permanent magnet eliminates the 
rotating armature and the mechanical wear normally associated with brushes. 
These motors fall in the class of Permanent Magnet Motors and enjoy certain 
advantages over wound-field types such as: 

"...Linear torque-speed characteristics , high stall (accelerating) 
torque, no need for electric power to generate the magnetic flux and 
a smaller frame and lighter motor for a given output power" [1]. 

Additionally, the brushless dc motor is characterized by: 

"...controllability over a wide range of speeds, capable of rapid 
acceleration and deceleration, convenient control of shaft speed and 
position, no mechanical wear problem due to commutation and better heat 
dissipation arrangement" [1]. 

The fundamental requirement of an electro-mechanical actuator control 
system is to provide torque to an output shaft, sense the position of the 
shaft and adjust the torque to balance the load when the desired position is 
reached. This must be accomplished with a minimum of frictional resistance 
and delays associated with the inertia of the mechanical components. Effects 
of viscous, static and coulomb friction, together with the torque required to 
accelerate the mechanical components of the system, lead to a reduction in 
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torque available at the output shaft and an associated reduction in system 
performance. 

One approach to the analysis of the electro-mechanical actuator system 
has been to divide the system into two sequential problem areas. The first 
deals with the dynamic analysis of the brushless dc motor and development of 
the transfer function necessary to duplicate actual steady-state and transient 
performance. The second area deals with modeling the mechanical system 
elements, taking as input the dc motor shaft angular acceleration predicted by 
the motor analysis. The mechanical system must be modeled considering the 
effects of friction and inertia and translating the rotational motion of the 
brushless dc motor shaft to the output shaft of the actuator for application 
to missile maneuvering control. The results obtained from a study that has 
placed primary emphasis upon the latter problem area - the modeling of the 
mechanical drive-train leading to the fin shaft was presented in Ref. 3 an d 4. 

This report deals with the first problem area and presents an overview of 
the work done by MacMillan [H]. In Ref. MacMillan developed a model for 
the output circuitry of a transistorized power conditioner that provided the 
required motor commutation. The model developed in Ref. 4 used an ideal (zero 
impedance) power supply. This report extends the model to include internal 
power supply resistance and also identifies areas for further improvement in 
the system modeling. 

The next section of this report presents a brief description of the 
system followed by a development of the model and an analysis of the simulated 
motor performance. Results of the simulated system are then compared to the 
measured response from a typical commercially available motor and recommen- 
dations for further improvements in the model are outlined. 
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SYSTEM DESCRIPTION 



GENERAL 

The system is viewed as a position control device to maintain an output 
angle under an applied hinge moment due to aerodynamic forces on a fin or 
aileron. The motor is a permanent magnet dc motor with feedback in the form 
of back emf proportional to the angular velocity of the motor. The block 
diagram of the dc motor and power conditioner is shown in Fig. 1 [4], The 
mechanical actuator and drive train, as currently envisioned, introduce 
various inertial and damping loads together with an aerodynamic force and its 
associated fin hinge moment that must be overcome to produce output motion. 
An operational block diagram of the load torque is shown in Fig. 2 where the 
hinge moment and motor shaft angular acceleration are viewed as inputs to the 
drive train [2], Figure 3 is a schematic of the drive train which, as 
presently constituted, includes the motor shaft (leadscrew), ball screw 
assembly, and the crank which is keyed to the output (fin) shaft. Inertial 
loads are considered individually within three major subdivisions of the 
actuator; the output shaft to crank, crank to ball screw, and ball screw to 
leadscrew. 

A detailed analysis and development of the mechanical model with the 
associated assumptions are contained in Reference 2. As stated in the 
introduction, this report reviews the development of the motor model reported 
in Reference 4. For convenience in the development of the model, MacMillan 
considered the load seen by the motor to be caused by the inertia and kinetic 
(viscous) friction of the motor shaft as described by Thomas in Reference 5. 

The schematic diagram of the power supply, power conditioner and motoi — 
load is shown in Figure 4 where it is noted that the power conditioner logic 
and driver circuits are modeled as ideal on-off devices with zero time delay. 
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Figure \ Motor and Power Conditioner 

Simulation Diagram 







Figure 2 Load Torque Block Diagram 
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Figure 4 



Schematic Diagram of Power Supply, Power Conditioner 
and Motor with Inertia and Viscous Friction Load. 
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Additionally, MacMillan considered the power supply to be ideal (zero interval 
resistance) to further reduce the complexity and simplify the development of 
the power conditioner model. However, in this report, the model used for the 
power supply includes the internal resistance, R s . The effect of including R s 
in the requirements and performance of the system model is presented in the 
section on Analysis of Simulated Motor Performance. For additional informa- 
tion on the development of the D.C. motor model refer to Reference 4. 

DEVELOPMENT OF SYSTEM MODEL 

The basic simulation by Thomas used a single power supply and super- 
imposed the phase currents to produce the motor torque. Motor drive was then 
realized by multiplying averaged armature current by a torque factor. 

In the model developed by MacMillan, the supply is considered to be a 
split supply of equal voltages. Armature current is assigned a positive sign 
if it flows in the positive direction (i.e. into the motor). The use of a 
split power supply allowed MacMillan to apply circuit reduction techniques. 
Because of the symmetry that resulted from the split power supply approach, 
the complex 3~phase bridge-type circuit simplified into a two-window network. 
The development of the reduced circuit begins with the definitions of loop 
currents as shown in Figure 5 where the power supply voltage of 2 V is 
represented as V+ and V-. The simulated model equivalent circuit is shown in 
Fig. 6 with node N defined as the mid-point of the power supply and no'de 0 to 
be the mid-point of the motor windings and neither of these nodes are 
considered to be at ground potential (the ground is identified as the 
Reference node). Note that the power supply, VIN = VIF - VIB = 2V and that 
all the network variables are defined using the CSMP model variable 
definitions [6]. Given the assumption that the network is balanced, that is, 
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Figure 5 Motor and Commutating Transistors 
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all transistors, diodes and field windings are approximated as being alike in 
characteristics, it is possible to apply the circuit reduction technique of 
Thevenin and obtain the simple 2-window network shown in Figure 7 [4]. In 
Figure 7, the Thevenin equivalent voltages and currents are defined in terms 
of the CSMP model variables. 

The back emf voltages of each phase (VEMFA, VEMFB and VEMFC) are 
developed and presented in Reference 4 in which it is shown that these 
voltages are summed 2 at a time in proper sequence to compute the loop 
currents and resulting phase torque. The total developed torque is then the 
sum of the 2-active winding phase torques over the proper 60 degree of 
mechanical angle. Figure 8 shows the back emf voltages generated across 2 
windings taken 2 at a time. The 60 degree of mechanical angle that produces 
torque is shown by the logic unit output levels of the position transducer 
RPSA , RPSB and RPSC in Figure 8. 

Reference 4 presents a detailed development of the power conditioner 
model that includes the assumptions for the switching transistor dynamics used 
in the power conditioner model as well as the development of the harmonic 
air-gap flux used in the motor model. 

ANALYSIS OF SIMULATED MOTOR PERFORMANCE 

The proper evaluation of a system model requires, quite naturally, a 
comparison of the actual system response for the same type of input function 
used for the model. For the brushless D.C. motor, it is assumed that 
measurements can be obtained for the following voltages and currents: 

- Power supply terminal voltage; V sg 

- Voltage across two of the three windings; V a ^, V^ c , and v ca 
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Figure 5 Split Power Supply Bridge-Type Circuit Model 
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Figure 6 Two-Window Equivalent Circuit of the Balanced Bridge Network 
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Fiqure 8 Back EMF and Position Sensor Output with 
Driven at 1200 RPM (CCW) 
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- Voltage from one winding terminal to ground; V a . V b and V c 

- Power supply current; I m 

- Phase current; I ma , I mb> I mc 

Before taking data from the simulated model, it is necessary to set all 
constants and parameters to the nominal values given by the supplier of the 
motor. Care must be taken to properly interpret the values given by the 
manufacturer . 

In particular, the value of the Back EMF coefficient, k b , needs special 
attention. If the motor data sheet indicates that this coefficient is 
obtained by measurement of peak to peak voltage across two windings, then in 
the model that uses distributed sinusoidal air-gap flux, an adjustment factor, 
k adj is required to insure that the Back EMF values are in agreement. It 
follows that the torque coefficient, k b , must also be adjusted by this same 
scale factor [7]. The value of k ad j can be computed from the measured no-load 
data for the motor current, speed and applied voltage by using the average 
characteristics balance equation, 

*m = ( v sg ~ k b w nl^^a + ^q) 
where k b = k bm k adj 

k bm = manufacturer supplied Back EMF constant. 

Thus the adjustment factor becomes 

k adj = ^ v sg ~ i-m^a + Rq)]/(k bm w n ]_) 

For example, given a no-load speed of 3060 RPM with a terminal voltage of 
30 volts and a no-load current of 0.30 amperes, k ad j = 0.825 (given a total 
series resistance of 1.47 ohms). The value of k ad j must of course be less 
than unity, otherwise the motor current would be negative in value which 
implies that generator rather than motor action is taking place. 
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Before an attempt is made to close the loop on the motor-load to form a 
positioning device, it is important to operate the system as a velocity device 
and obtain measurement of the motor average performance characteristics. 
Typical motor average performance curves are speed versus torque, motor 
current versus torque and output power versus torque from no-load torque to 
near stall torque conditions. Figure 9 shows these three curves for the model 
with torque given in oz-in units. Note in particular the straight line 
characteristics for motor speed and current and that the power output is of 
quadratic form. The condition of constant input voltage is used when 
gathering data for these curves, and it follows that if speed is a straight- 
line versus torque, then current also has a straight-line when plotted against 
torque. This can be shown as follows: 

w m = w nl - k wlT 
where w m = motor speed 

w nl = n O“load speed 

k w i = magnitude of slope of the speed vs torque curve 
T = load torque 

The current, using average values, can be written as 

*m = ( v sg “ k b w m^^a + ^q) 
where v S g = power supply terminal voltage 

k b = Back EMF coefficient 
R a = winding resistance 
Rq = transistor "ON" resistance 
then by substitution, 

i-m = t( v sg ~ k b w nl^ ( R a + R q)J + C (k^ky^ )/ (R a + Rq)]T 
where the first term on the right hand side of the equation represents the 
no-load current. 
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Figure 9 Average Performance Characteristics of the Simulated Model 
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The power output versus torque curve in Figure 9 reaches a peak value 
between no-load torque and stall torque. The analytical development for this 
quadratic curve is presented in Appendix A where it is shown that the 
condition for peak power output, P 0 , is that the torque must be of value 

T = (2k^w n ^ - v S g ) / ( 2k|- ) k w ^ ) 
to produce maximum power output 

P 0 = k t v sg 2/ ^ R s k b^ 

Since the ratio k^/k^ is a constant, increasing peak power output depends 
upon decreasing the power supply interval resistance, R s , or increasing the 
power supply voltage, v sg as one would logically expect. It must be 
understood that the assumption used in the above analysis was that the air- 
gap flux was a constant average value. For the actual motor the air-gap flux 
is distributed sinusoidally and current in the windings produces a field that 
may cause distortion in the air-gap flux. At large current values that occur 
at and above peak power output, the field distortion may result in an increase 
or a reduction of output power for a given load torque. 

Another motor coefficient to consider carefully is the no-load viscous 
friction constant, B m . The value of B m is often not included in the motor 
data sheet and the probable reason is that its value will depend somewhat on 
the manner in which the motor is attached to the system load. B m can be 
calculated from given no-load data as follows: 

B m = ( k t k adj i m) /w nl 
where i m is the no-load current. 

RESULTS AND CONCLUSIONS 

The preliminary results indicate that the balanced bridge circuit 
approach used in the development of the Brushless DC motor model produces good 
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agreement with measured motor characteristics as indicated below. Further 
evaluations of the model as well as improvements and additions to the model 
will be conducted in the near future using data currently being gathered from 
a prototype of a fin positioning actuator. 

A comparison of model and motor produced Back EMF voltage is shown in 
Figure 10. The upper curve is the simulated motor waveform and the lower 
curve is the measured Back EMF. The model adjustment factor, k ac ij was set to 
0.63 value by trial and error until peak to peak voltages were in agreement. 
Another verification of the model is obtained by comparison of Figure 11 and 
Figure 8. Both figures show waveforms of Back EMF across 2 windings taken 2 
at a time and also show the timing waveforms for the Position Sensor Devices 
( RPSA , RPSB and RPSC) for counter clockwise rotation. The waveforms agree in 
both phasing and in form. 

The steady-state performance curves of a typical motor are given in 
Figure 12 where the motor load torque is given in lb in units. These curves 
were produced with a constant terminal voltage of 30 volts and agree in form 
with the curves of the model as given in Figure 9. Peak power output occurs 
at load torque of approximately 95 oz-in for both motor and model. 

Additional validation of the model is shown in Figure 13 where the 
current in Phase C for the model (upper curve) is in close agreement both in 
form and in phase with the same current for the motor (lower curve). 
Additional model results are shown in Appendix B where typical input voltage 
step response waveforms are presented. 



FURTHER RESEARCH 

Improvements in the model are required with regard to the triggering on 
of the diodes used to protect the transistors from excessive reverse currents 
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Figure 10 Back EMF Voltage: Top Curve-Model Output, 

Bottom Curve-Typical Motor Output 
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Figure 11 Back EMF and Position Sensor Output Waveform for a Typical Motor 
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daring the switch off time for the transistors. MacMillan in Ref. 4 was able 
to provide this action for the condition that the power supply interval 
resistance, R s is zero value. When R s is not zero, additional programming is 
required to perform the required trigger action. 

The Pulse Width Modulation (PWM) used by Askinas in Ref. 8 must be added 
to the present model, however Askinas used a separate additional transistor to 
perform PWM. The Power Conditioner to be used in the Test Stand for 
measurement of actuator performance accomplishes PWM by using the lower set of 
the commutating power transistors. Therefore some modification of Askinas’ 
CSMP coding may be required to accurately perform the PWM speed control. 

Extensions of the model will include a position sensing device and a 
tachometer as well as the nonlinear loading as developed by Franklin in Ref 9. 
Test data results will be compared to model in both transient (step) response 
and in frequency (Bode) response and the results will be reported in a Summary 
Report for FY1986. 
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APPENDIX A 



Power Output of a Brushless DC Motor 

The mechanical power developed by a Brushless DC Motor is given by the 
equation, 

P = (7 .062x10“3 ) Tw m watts 

where T=Tl+Tm, the sum of the load torque, T L and the motor restraining 
torque, in oz-in units and w m is the motor speed in rad/sec. units. At no 

load Tl = 0 and at stall torque w m = 0, thus power reaches a peak value in the 
region between no load and stall torque conditions. 

The evaluation of steady-state peak power can be obtained analytically in 
terms of the load torque by substitution as shown below. 

Given that T=Tl that is T[ j >>T[ U ] and T = k t i m with w M =w f j[_ i where k b is the 
motor torque coefficient and k w ^ is the slope of the speed-torque curve. 

Then the power output is 
P 0 = kt^m ( w ~ ^wl^L^ 

and since steady-state speed conditions are assumed, the current is given by 

v - k w w 
. _ sg b m 

m " RA + Rq 



where v sg is the voltage at the terminals of the power supply, k b is the Back 
EMF coefficient, R a is the armature resistance of 2 phases in series and Rq is 
the sum of 2 power transistors forward conduction resistance. By 
substitution, 



k. v k. k, ; 

P = — - — f w - |< T 1 - — - — — fw - k T 1 
0 R a + R q nl wlV R^+Rq [ nl wl L J 



and to find the conditions for peak power output, it is necessary to compute 



A1 



the rate of change and set it equal to zero. Thus, 



dT L = ' R A + Rq 





Solving for load torque, TL yields 

Tl = w nl^ k wl ~ v sg/ ^ 2k g k wl ) 

T l = (2k b w nl - v sg )/(2k b k wl ) 

and substitution into the power output equation produces the conditions for 
peak power output as follows: 

F o = k t v sg / ( R a +R qHwnl _k wl(2k b w nl -v sg )/(2k b k wl )] 

~ k t k b // ( R a +R q) [ w nl~ k wl ( 2k b w n ^-v sg ) /2k b k w ^ ) ] 2 
After expansion and cancellation of like terms, the peak power output reduces 
to 



As a numerical example consider the following values of a typical 
Brushless DC Motor operating at a steady-state speed of 3060 RPM (320.4 
rad/sec) with a constant terminal voltage of 30 volts. 
k b = 15.9 oz-in/ampere 
k b = 0.112 volts/rad per sec 
k adj = 0.825 

R a = 1.37 ohms, Rq = 0.10 ohms 
k wl = 1 *67 rad per sec/oz-in 
Thus for this example motor, 



P Q m k b v sg /[4k b (R a +Rq) ] 



P Q = (15.9)(0.825)(30) 2 (7.062x10“3)/[4(.112)(0.825)(1 .47) 



P Q = 153.5 watts 



and occurs at load torque 



Tl = [ ( 320 . 4) /( 1 .67) ~ (30.)/(2)(.112)(.825)(1.67)] 
T l = 94.65 oz-in = 5.92 lb-in 



A2 



The curve of power output versus load torque for this example motor is 
shown in Figure 9. It must be noted that the above analysis is valid for the 
assumption that motor current produced fields do not have a measurable effect 
on the magnitude and phase of the flux generated by the rotating permanent 
magnets and also that a constant average flux exists in the air-gap between 
rotor and stator. 
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APPENDIX B 



CSMP Listing 

The system described by the schematic diagram in Figure ^ is programmed 
in the CSMP language as shown below. The computer used was an IBM 370 located 
in the W. R. Church Computer Center at the Naval Postgraduate School. 
Following the program listing is a User's Guide for the recommended procedure 
in setting the constants and parameters for validation with measured motor 
performance. Also included are typical output waveforms for a step input of 
supply voltage. 
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FILE: BMIRA 



CSMP 



A 1 



//JRBTB001 JOB (C 169, C 3 23) , ’MACMIILA-IB ' ,<-LABS = G,MSGLEVEL = ( 1 , 1) 

//♦MAIN CSG=NPGvM1.C1fc9f,L3NLS=SSS 

//♦FOR MAI EE ,DDNAME=,DESI=I C c A L 

//♦FORMAT EE , DD N AM E =S YS V ECT R , L E S T = 1C C AL 

//CSMPSTEP EXEC CSMP VD V , POE DI E S = K Y C 1 E S 

//X.SYSIN CE ♦ 



♦ PATCHES E BOB IOC INSTALLED 

♦ VERSION ELEVEN -- SIXTH REVISION I NCCE ICR ATI NG IHE3V5 LIU. 

♦ THIS VERSION INCORPORATES EEVEESING COMMUTATION AT EEVTIM. 

♦ CURRENTS ABE NCI SUPERPOSED. ELUX * CURRENT IS COMPLIED FOR EACH EEASE 

♦ AND THE RESULTING TORQUES ARE SUMMED. 

♦ THE PURPOSE OF VERSION FIVE IS TO TREAT THE EIUX AS VAFYING 

♦ ACCORDING TO THE SUM OE A FUNDAMENTAL SINUSOID 

♦ AND ITS FIFTH HARMONIC AS EXPLAINED IN CHAPTER FIVE. (THOMAS) 

♦ THE TOTAL FLUX IS AFPECXI BATED AS THE AIGEEFAIC SUM CF THE EIUX 

♦ DEVELOPED IN TWO W I N DINGS AT A TIME IC CALCULATE THE 1C5£UE 

♦ COEFFICIENT AND LACK EMF CCEEEIC1ENT. 

♦ THIS MODEL SIMULATES MAGNETIC FLUX AS A SINUSOID AND A EIETB 

♦ HARMONIC AND SIMULATES DIODE COMMUTATION AS WEIL AS IHE SNITCHING 

♦ LOGIC AND TRAN SIS I OR DYNAMICS CE A BRUSHLESS DC MCTCR POWER 

♦ CONDITIONER. THE WINDINGS ARE NCT TREATED INDEPENDENTLY 

♦ AND THEIR TRANSIENT ELECTRICAL I NT h R F E I A 1 IC N IS SIMULATED, 

♦ THE CONTRIBUTIONS 1C DEVELOPED TCF£UE ARE 

♦ TREATED AS SUP ERPCSAEI E . 

♦ 

* THE PROGRAM HAS BEEN MODIFIED 10 INCLUDE POWER SUPPLY E ESI ST ANC £ , E S . 

♦ WITH THE ADDITION CF B3, THE 1PIGGEE1NG VOLTAGE ECR THE TRANSISTOR 

* PEG1EC1ING DIODES HAS ££EN CHANGED TC INSUEE P EC F ED CEEFA1ICN. 

♦ 

♦ 

* 

INI1IAI 

♦ 

CONSIANT KTK = 15.8b, EM = 0.01485, El = 0.00, JL = 0.0, N = 1.0, ... 

JM = 0.C01. REM = 0.1120, PI = 5.14159265, K K 5 = 3.59C,KADJ = .6*, 

THA DV = 0.5, ?S= 8.53, IC1MC = C.C1 

* 

PARAMETER LA - .C008, FA = 0.6S5C, TII= C.C 
PARAMETER VSAT = .4, RSA1= .05, FOCI = 1.CE+4 

PARAMETER ITEMS = 1.0E-6, REVIIM = C.10 

PARAMETER VC1D=0.,VC2C=C.,VD3D=0..VC4C-0.,VD5C=C.,VC6E=C. 

PAKAMATEB VSGPEE = 3C.C, KINT = 1CC0C.0, KRAKP = 18250. C 

+ 

* 

♦ KTM — MEG SUPPLIED MOTOR TORQUE CONSTANT (CZ-IN/AMP) 

♦ KT — PHASE ICBCUE CONSTANT (CZ-lN/AflE) 

♦ KBM -- MEG SUPPLIED MOTOR HACK EMF CONSTANI ( VC IT/ E A D/S EC ) 

♦ KE — PHASE BACK DBF CONSTANT (VCIT/RAC/S) 

♦ KADJ -- ADJUSTMENT EACTOR FOR MEG SUPFIIEL MCTCF CCNS1ANIS 

♦ RA — PHASE RESISTANCE OE THE MCTCF (OHM) 

♦ LA — PHASE INDUCiANCE Cr THE MCTCR (HENRIES) 

♦ 3M -- VISCOUS FFICTICN COEFFICIENT CF THE MCTCR (C Z - 1 N/ R A E/S ) 

♦ BL — VISCOUS FRICTION COEFFICIENT CF THE DC A D (C Z- 1 N/ E A E/ S) 

♦ 3LP — VISCOUS FRICTION COEFFICIENT CF LOAD THRU REDUCTION GEARS 

♦ E — TOTAL VISCOUS FRICTION CF THE MOTOR SYSTEM (SEE BELOW, NC - S"C Rl) 

♦ JM -- .iNEETIA CF THE HCTC* (CZ-1N/S-S) 

♦ JL — INERTIA CF THE LOAD " " » 

♦ JLP — INERTIA OF THE LOAD THFU REDUCTION GEARS 

♦ J — TCIAL INERTIA CE THE MCTCF SYSTEM (SEE BELOW, NC-SCRT SECTION) 

♦ A2 = J/E — THE MECHANICAL TIME CONSTANT CF THE MCTCR (SEC) 

♦ ABTAU = LA/ (RA+EEQAE) -- IHE EIECIPICAI TIME CONSTANT CF THE MOTOR 

♦ ANC Lb I V c TRANSISTORS FOR CURRENT PATH A-E 

♦ BCTAU = LA/ (RA+REQDC) -- THE EiECTblcAI TIME CONSTANT CF THE MOTOR 

♦ ANC DRIVE IRANSx STOPS FCF bUFREM PATH E-C 

♦ CATAU = L A/ (SA+REQCd) -- 'I li E ELECTRICAL TIME CONSTANT CE IHE MCTCR 

♦ AND DRIVE TRANSISTORS FCF CURRENT PATH C-A 

♦ V I F = 1/^ SUPPLY VOLTAGE, VIN, MEASURED FROM + TC MIC £CSIIICN,N 

♦ V IB = 1/2 SUPPLY VOLTAGE, VIN, MEASURED FRCM - IC MIC ECSIT1CN,N 

♦ HQ# -- EQUIVALENT CIRCUIT RESISTANCE CF Cl*) 

♦ Q#EXP = E X PO j £ N 1. 1 A I MODEL FOR TFAaSIS'UCN BETWEEN SATURATION AND CUTCFE 
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C SME 



A 1 



ESAT = EQUIVALENT CIRCUIT PE3I3TANCF. CE 1EANSI3ICR SA1UFATICN 
RCU1 = EQUIVALENT CIRCUIT RESISTANCE CE 1BAASISICR CUICEE 
TIIME — SWITCHING TIME CF TFANSIS1CE 

1HHSI —THIS IS IHDEG WHICH IS SET TC 2 £ F C AT STA5I, BESET 1C 2EFC 
AI 3o0 DEG AND CAN AISC EE USED AS SIMUIATICN SICE ANGLE. 
THADV -- THETA ADVANCE FOR COMMUTATION (CCW IS POSITIVE) 

THCCN — COMMUTATION CONTROL A N G L £ = 1 H F S X ♦ 1 H A C V IN DEGREES 
fiS INTERNAL RESISTANCE Cf 1HE SUPPLY VC11AGE 



DSCST 

KT=K1M*KADJ 



KB=KBM*KA£J 
BLP = EL/ (N**2) 

JLP = J L/ ( N * * 2 ) 

J = JM ♦ JLP 

E = B il ♦ EIP 

A 1 = LA / BA 

A 2 = J / E 

A3 = LA/ (BA +RSAI) 

RS1=E5/2. 0 

ES2=BS/2. 0 

VINIC = 15.0/KINT 



21HCE SI1EF 
(NAMIC 

VIF= VINHAE 

VIB =- V I N H A F 

VIN = VIE - VIB 

VSGDEL = (VSGREF - VSG)/2.0 

VSGEER = DEADS? ( -C.5, <-3.5, VSGDEL) 

VININT = INTGBi (VINIC, VSGEuR) 

VINHFG = K IN I * VININT 
VINHAF = VI.NEF1 - VINEF2 
VRAMP1 * EAMP(C.C) 

VRAMP2 = E AMP (3.00 1) 

VINHF1 = KRASF*VFAMP 1 
VINHF2 = K R A M f + V F A M P 2 



INTRODUCE UNITE TRANSITION TIME FCH IBANSISIOE SWITCH 
BY INCCFPC RATING EXPONENTIAL RISE AND DECAY IN1C SW1--SW6 
Q # EX P PROVIDES A REALISTIC EXPONENTIAL IFANSITICN EE1EEEN t.0TGFF 
AND SATURATION -HIGH INCLUDES SATURATION DELAY WHEN PDI 1HFCUGH 
THE LIMITER 



wlEXE = SEAL PL (.75, 11 I ME, 
Q2EXF = RE ALPL (0.75, I LIME 
Q3EXF = FEALPL (3.75, TIIME 
Q4EXP = iiEALPL <0. 75, TIIME 
Q5EXP = S £ A L ? L (0. 75, TIIME 
QoEXP = REAL?! (0.75, 1 X 1 3 Z 
:LE3E Q1EXP = J.1, Q2fcX?=C.1, 
3SEER Q 1 £ X P = 0 . 1 , £eEXP = 0.1, 
RQ 1 = LIMIT (ESAT, ROUT, 2. 
Rw2 = LIMIT ( f 5 A 1 , LCD I, 2. 
RQ3 = LIMIT (ESAT, ECU'!, 2. 
HQ4 = LIMIT (ESAT, ECUI, 2. 
RC5 = LIMIT (ESAT, RCU1, z. 
Ego = LIMIT (ESAT, PCUT, 2. 



3.75-SW1) 

, 0.75-SWz) 

9 



3 



0. 75-S 
0. 7 5 - S W 4 j 
, 0.75-SK5 
, 0.75-SW6 
Q3FXP=C. 1 
Q 3 E X t = C . 1 ; 
*ECU1*C1 EXP] 
* E C 0 I 4 v 2 E < F 
+ECUT*C2EX? 
+FCJT+C4EX? 

+ ECUT*\,5FXP 
♦ECOT+CcEXS 



*'4£Xr = 0. 1 
C4EXP=C. 1 



£5£XP = C. 1 
£ 5 5 X P = 0 . 1 



£6EX 
3 c E X 



LINEAR CICDE MODEL 
DIODE T0RN-CN VCLIAGES 



VD ID = 


6CC0.0 - VA.IND «• 


VCINE 




VD6D = 


bCGC.O + VCINE - 


V £ I N C 




VD3D = 


fcCCC.C - V ii I aC <- 


V A I N C 




VD2D = 


6COO.O + VAINE - 


V C I N C 




VD 5D = 


C3CC.0 - VCINE ♦ 


V E- 1 N D 




VD4D = 


cCOO.O + V E I N C - 


V A I N C 




R D 1 = 


EC NSW (VD1D, RSA1, 


SCUT , 


ECUI 


ED2 = 


F C N S W ( V D ^ C , ESmT, 


E C 'J X , 


ECL : 


RD3 = 


EC NSW (VB3E, R S A 1 , 


SCOT, 


ECUI 


8C4 = 


ECNSW (VD4D, RSAT, 


ii C O 1 , 


cC U 1 


F. D5 = 


ECNSW ( V D 5 C , ESAT, 


SCUT , 


FCU i 
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A 1 



R Do = FCNSW ( V C o C , FEAT, RCUT, rClI) 

* 1 HE IHEVSNIN E Q DIVALENT C I ruj (J I 'I METHOD IS USED WITH IRE SUPPLY 

* VOLTAGE MID-POINT , N, AS REFERENCE AND 1HE LCCP MEIriCC CF ANALYSIS 

* APPLIED IC COMPUTE THE PHASE CIFFFMS. 

* 



* THEVENIN ECUIVAI EM VOLTAGES 

VANTH = VN1*FEQAS2/ (f> EQ AS 1 -» F EQ AS..) 
VBNTH = VNl*E£QBS2/ REQBS1fREQES2) 
VC NTH = VN1*FEQCS2/ F EQCS 1 +F £ Q C S 2 



+ VN2*i»ECAS 1/ (REQAS UFEQAS2) 
+ VN2+FEQES 1/ (RFQtS 1 + FEQES2) 
+ VN2 + FEQCS 1/ (BEQCS 1 + REQCS2) 



TRANSISICr/DICCE EQUIVALENT 



REQA 1 
REQA2 
REQ E 1 
REQB2 
EEQC 1 
REQC2 



- B C 1 

- R Q 2 
= FQ3 
= RQ4 

- HQ 5 
= RQ 6 



SCI/ (FQ 1+RC1; 
P C 2/ ( F Q 2 + a D 2 
RCJ/ (FQE + ^Dj' 
RE 4/ (FS.4 + F.C4 
EC 5/ (FQ5 + RE5 
RE6/ ( f v 6 ♦ R D 6 ] 



RESISTANCE 



* COMPUTATION CF IHEVENIN EQUIVALENT RESISTANCE 

* 

REQ AS 1 = REQA1 + RSI 
BEQAS2 = REQA2 + F S 2 
REQBS1 = FE*E1 + RSI 
REQBS2 = FEQE2 ♦ F S 2 
REQCS1 = REQC1 ♦ RSI 
REQCS2 = FEQC2 + R S 2 

* IHEVENIN EQUIVALENT RESISTANCE 

REQA = FEQAS1 * FEQAS2/ (REQAS1+FEQAS2) 

REQB = EEQES1 * F E w BS 2/ ( F E w BS 1 + F EQ LS 2) 

REQC = REQCS1 * R EQCS 2/ (R E QC S 1 ♦ R E Q 0 S 2 ) 

* 

* TOTAL LOOP RESISTANCE 

REQAE = 2 . *R A + F.EQA ♦ EEQE 

REQBC - 2.+HA ♦ REQt + REQC 

REQC A = 2 . * P A + IEQC * REQA 

* 

* LOOP CURRENT- MAXIMUM VALUE 
IAB - (VANTH - VENTF) /REQAt 

IAB A = ♦ ( V I N + V E M Ft: - VF.MIA)/(FS ♦ 2.*3A + FEQa2 
IBC = (VENTH - VC NT fi ) / REQBC 

IBCA = -(VI N + VEXES - VEMFC)/(FS + 2.*FA + REQC1 
ICA = ( VCNT H - V A NT E ) / 5 EQC A 

ICA A = - ( V 1 N + VEMFC - VZMEA)/(RS + 2.+RA ♦ FEQA1 

* 

* TIME CONST ANTS 

A El A U = 2.* LA/REQAP 
BCTAU = 2 . * LA/EEQBC 
CATAU = * LA/FEQCA 

* 

* LOOP CUFF ENTS 

IMAB = 5EALPL (0. 0 , AET AU , IAL) 

1M3C = REAL PL (C.G # EcIAU,IcCj 
IlMCA = FEALPL <0. 0, CATAU, ICA) 

* NET LEG (PHASE) CURRENTS 
IMA = I K A fc - I MCA 

IMB = I M Be - IMAB 
IMC = I MCA - IMEC 
IMCD = DERIV (3. 0,IMC) 

IMCCR = TCIMc+IMCC 

IMCAC = hEALPL (0. 0 , TC I MC , I M.CDF ) 

IMF = REAIPL (0.0,0. CC5.IM) 

VSGF = FEALPL (0.0,0.005, VSG) 

* 

* INDUCTOR INDUCED VOLTAGES 
CIADT = DFuIV (Q.C, IMA) 

DIBDT = DERIV (C.0, I E ) 

DICDT = DERIV (C. 0, IMC 
VA1NC=LA*I I AD I 
Vt3INC=LA*CIt:DI 
VCINC=LA+DIOD I 

* PHASE OEM'S L A « VCLIAGF 



+ F E Q A 1 ) 
♦ BEQE2) 
+ FEQC2) 
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L E : E a *1 B A 



CSHP 



A 1 



V AOL V - 
VBOL V ; 
VCOL V : 
PHASE 

V AO = 

V BO = 
VCO = 



im a* ( r a 

I K E * ( B A 
IKC* (BA 
VCLlaOE 
ISA* (HA 
IHB* (BA 
IMC* (BA 



♦LA*CI ACT+VEMFA 
♦ LA*EIBCT + VEf'. F£ 
+LA+EICEI+VFMFC 



COMPUTATION OF PCWEB SUPPLY TERMINAL VCLIAGES (GND BEE) 
WHERE VIS IS SUPPLY VOLTAGE A HE IH IS SUEELY CUHEES1 
VSG = VIN - IB+BS 

VAG = VSG * EEGA2/ (FECA1 + R EQ A2) 

VEG - VSG * EECB2/(BECE1 ♦ FECE2) 

VCG = VSG * EECC2/(EECC1 ♦ R EQ C 2) 

VAC = VAG - VCG 

VOLTAGES ACROSS TBANSISICHS ANC EICCES 

VQD1 = VSG - VAG 

VQD2 = VAG 

VQD3 = VSG - V EG 

V£D4 = VEG 

VQD5 = VSG - VCG 

VQD6 - VCG 



P V AG 


= 


LIMIT 


(-60. , 


6 C . # 


VAG 


P VBG 




LIMII 


-60., 


ec.. 


VEG 


P VCG 


X 


LIMIT 


(-60. , 


6 C . , 


VCG 


P VAC 


X 


EVAG- 


P VCG 






P V B A 




EVEG- 


P V AG 






PVCB 


X 


PVCG- 


P V EG 







Id, 

IQ2 
IQ3 
IQ4 
IQ5 
IC6 
PC 1 = 
FQ2 = 
PC 3 = 
PC4 = 
PC5 = 
pq6 = 
PQTOT 



XT 


FCNSJi 


(s-n. 


. 003 , 


. 002 , 


IMA 


X 


c C N S *> 


(S W 2 , 


.003, 


. 002, 


IMA 


X 


FCNSW 


(SH3, 


. CC3 , 


. 00 2, 


I ME 


X 


F C N S M 


(SW4, 


.003, 


. 002, 


I HE 


X 


FCNSW 


(SW5, 


« 0 C j , 


. 00 3, 


I KC 


= 


F C N S V 


( S ft b , 


. 003 , 


. 003 , 


IMC 



VEQ 1 
VCQ2 
VEC 3 
V D C 4 
VEQ 5 
VEQb 
= PC 1 



IQ I 
ICi 
IQ3 
IC4 
IQb 

I Se C 

PC2 



+ PQ3 + PQ4 ♦ EC 



PQfc 



THE PHASE ANGLES LCCK CCNEUSEC BUT 1 HE PHASE FELATICNSH1P IS RIGHT 

BEMFA = (3.*3IN(2.*lhEIA*(11*PI/t)) +.5b*SIN (1C.*?HETA*{1*Pl/6)) ) 

BEflFb = ( 3 . * S I N (2.*1HEIA+ (7 + PI/c) ) +.5y*SIN (10.+THE1A* (5*Pl/o)) ) 

EEHFC = i 3 . * S I N i2.*IHETA* ( 3 * P I / 6 ) ) +.59*SIN (1C. *T MET A* ( 5 * f I / b ) ) ) 

NCRMALI2E LEG £ M E 

VErtEA = BEMFA * K 0 / K K 3 * n M 

VEfIFE = EEMFE * KE/KK2 * HE 

VEMFC = B E MFC * KE/KK2 * WM 

VEHFAC = VEMEA-VEMFC 

V E ME E A = V2HFE-VEMFA 

VEHFC E = VEMFC-VEMFE 



TORQUE IC HAKE MCTCF TURN 

TA = IMA * KT * EEMEA/KK3 

TB = I HE * KT * EEMFE/KK3 

TC = IMC * KT * EEHEC/KK3 

TM = TA * Tfi + TC 

TMM = KI * IH 

TBH = EH * WH 

VKBWM = K E * W M 

TORQCE INTEGRAL VALUE 

THIN 1 = 1 NTGRL (0. 0,IM) 

TL = TLL+SIEP ( 1 . C £ - 3 ) 

TN2 = TNI * (1.0/E) 

W M = REALEL (C.0,A2,1N2) 

K MINT = INTGrL ( 0 . 0 , * M ) 
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* wm=32G.44P5 

WBRPM = WK * (20. /FI) 

WMRPMF = XMRFM/N 
THETA = IMGBi (O.O.hM) 

IHDEG = THETA * (IfiC.C/PI) 

* 

* COMMUTATION ADVANCE HEBE 
THCON = THRSI+THAO V 

* 

PWRML = Ail * T H * . 3C 7 3625 
PWRM = EM * .0070675 * WM + + 2 
PWR = fcL * . 3 0 7 0 c 7 5 * *M + *2 
IL = V IN/ (2*RCUT) 

* THIS PROCEDURE FRCVIEES A SIMPLE MECHANISM EOF REVERSING THE MOTCB'S 

* DIRECTION //HUE KEEEING THE LOAD AS AN OPPOSING TCECUE. 

PROCEDURE TN 1 = FWEEhC (WK,TM,TI) 

If (HB.LT.G.O) GC TC 1 
TNI = IB - T I 
GC TC 5 

1 TNI = 1M + TL 
5 CONTINUE 
ENDPRCCEEU RE 

* THIS PROCEDURE RESETS THE VAEIAELE THRST TC C AFTER EVERY 360 DEGREES 

* OF MECHANICAL ROTATION. THIS IS FUNDAMENTAL TC THE SIMULATION OF ALL 

* SWITCHING AND POSITION SENSING ACTION. 

PROCEDURE THRST = RESET (IFCEG) 

IHEST = A MOD ( I F C E u , 2oJ.) 

IF (THRST .11. 0.0) TORSI = THRST ♦ 260. 

ENDERCCEDURE 



* THIS PROCEDURE V, AS TAKEN FROM IhE NEXI PROCEDURE 

* SENSORS SHOULD GENERATE A SINGLE VALUED JUNCTION 

* REGARDLESS OF DIRECTION OF nCTATIGN 



PROCI 


•DURE SE 1 ,Sf^.SE2 = BALL (THCON) 








IF (THCON 


. GE. 130.) GC 


1C 4 5 










GC 


TO 46 














45 THCON = 


IHCCN - 180. 












46 CONTINUE 














IF I 


[T HOC N 


.GE. 0..AND. 


IHCCN 


•LT. jC.; 


| GC 


1C 


10 


IF | 


THCON 


• G £ • 3D**ii!sL/» 


THCCN 


. LT. 6 0.; 


> GC 


1C 


1 1 


IF | 


THCON 


• (j sL m 60* 


IHCCN 


.IT. 50. 


I GC 


1C 


12 


IF { 


THCON 


• G 2 . 9 C . .AND. 


IHCCN 


.LI. 1*0. 


1 G C 


1C 


1 3 


IF l 


THCON 


1^0.. A NC. 


IHCCN 


. L T . 1 5 U . 


1 GC 


1C 


1 4 


IF 1 


THCON 


.GE. 153. .AND. 


i H C C M 


.11.180.] 


\ G C 


1C 


15 


10 i 


;ei = 1 


• 















SEP = 


0 




SE3 = 


1 




GOTO 


PC 


11 


SE 1 = 


1 




SEP = 


c 




SE2 = 


u 




GOTO 


PC 


12 


SE 1 = 


1 




SEP = 


1 




SE3 = 


0 




GCTC 


PO 


13 


S E 1 = 


c 




SEP = 


1 




S E 2 = 


0 




GOTO 


PO 


14 


S E 1 = 


\J 




SEP = 


1 




SE2 = 


1 




GOTO 


<'■' •> 


15 


S E 1 = 


0 




S E * = 


0 




SE 2 = 


1 


20 


CC NT I 


MU 



BECAUSE THE HAIL 
Of POSITION 
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DFBCCECURE 

THIS PROCEDURE SIMULATES GENERAL COMMUTATION, 

DETERMINES IHE ALGEE5AIC SUM CF 

THE VAFIAELE FLUX CCMECNEN1S FCR USE IN COMPUTING THE MCICF'S 
APPROXIMATE SACK EM F (BEMFT). IT ALSC SUE16AC1S THE GENERATED 
VCITAGE F B CM THE PECEEE SUPPLY VCITAGE. 

SKI 1HEU S W 6 

SIMULATE POWER TRANSISTOR TRIGGERS EEING ENERGIZED CR SWITCHED OFF 
THCON IS THE VARIAELE THFCUGH WHICH THE SWITCHING ICGIC 
IS IMPLEMENTED. REVTIM IS 1 HE TIME AT WHICH CLOCKWISE COMMUTATION 
BEGINS . 



OCEDUEE 3Wl,3'W2,SW3,SW4,3W5,SW6,EEfF'I,VN1,VN2 =00 MM <11 f. E 
V1IM , THCCN , BEMFA, E EMF E , BEMFC , VEMF A, VEMFE , VEMFC) 

IF (TIME .GT. REVTIM ) GC 1C IOC 
THC0N1 = A MOD (THCCN, 1«0.) 

IF (T HCO N 1 .LT. C.C) IHCONl = THCCN1 4 180. 

"..AND. THCON 1 .LT. 3C.) CC 1C 

IHCONl .IT 
THCON 1 .LT 



IF (THCCN 1 
IF i THCON 1 
IF {THCCN 1 
IF (THCCN 1 
IF (THCCN 1 
IF (THCCN 1 



• GE. 0 . . A N D ■ 
,GE. 30. .AND 
.GE. bG . . A ND . 
.GE. 9 0 . . A N C , 
.GE . 120. . AND , 
.GE.15U..ANC, 



THCCN 1 .ET 



it TiH-Lfl I .bt.l~U..fl« 

CLOCKWISE COM MUTAIICN 



THCCN 1 
THCCN 1 



. LT 
.IT 



4 

30. 

60. 

SC. 

120 . 

150. 

180. 



GC IC 
GC 1C 



50 

5 1 



GC TC 53 



GC 

GC 



1 C 
TC 



54 

55 



0 CONTINUE 
IF (THCCN 1 .GE. 
IF (T HCC N 1 .GE. 
IF (THCC N 1 .GE. 
IF i T HCC N 1 .GE. 
IF (THCCN 1 .GE. 
IF THCCN 1 .GE. 
G SW 1 = C. 



0 . .AND, 

30.. AND, 

60. . AND. 

50. . AND, 
12C..ANC, 
150. .AND, 



SW2 = 


0. 


3 W 3 = 


C. 


SW4 = 


1. 


3W5 = 


1. 


SW 6 = 


0. 


BEMFI 


= HEMFC - E E M f n 


VN 1 -= 


V I F - VEMFC 


V N 2 = 


VIE - VEMFE 


I M = I M C 




GC IC 


60 


SW 1 = 


1. 


SW 2 = 


c . 


SW 3 = 


0. 


S W 4 = 


1. 


SW5 = 


c. 


S W 6 - 


0. 


BEMFI 


= EEMFA - EEMFB 


VN 1 = 


V I F - V E M f A 


V N 2 = 


VIE - VEMFE 


I M = I M A 




GC TC 


60 


SW 1 = 


1. 


SW2 = 


0. 


SW3 = 


0 . 


SW4 = 


0. 


SW 5 = 


c. 


SW 6 = 


1. 


BEMFT 


= BEMFA - EEKFC 


VN 1 = 


V I F - VEMFA 


V N 2 = 


VIB - VEMFC 


I M = I K A 




GC 1C 


cC 


SW 1 = 


G. 


SW 2 - 


0. 


SW 3 = 


1. 


S W 4 = 


0. 


SW 5 = 


c. 


S W 6 = 


1. 


BEMFT 


= EEMFE - EEMFC 



'IHCONl 
THCON 1 
THCON 1 
THCON 1 
THCCN 1 
THCON 1 



IT. 30.) 
L 1 . b 0 . ) 
IT. SC. 
LI. 120. 
11.150. 
LT .130. 



GC 


TC 




GC 


TC 


54 


GC 


TC 


C. C 


GC 


TO 


50 


GC 


TC 


5 1 


GC 


1C 


52 
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FILE: E .1 T & A CS.1E A I 



VN 1 = VIF - VEMEE 
VN2 = Vic - VE1FC 

i a = i h e 

GC 1C 60 

54 SW1 = 0. 

SW2 = 1. 

J H J la 

S W 4 = 0. 

SW 5 = C. 

3W6 = 0. 

BEKF'I = BEilfE - t E M E A 
VN 1 = VIF - VEMrE 
VN 2 - VIE - VEMFA 
IM=I ME 
GO 10 cO 

55 SW 1 = 0. 

SW 2 = 1. 

SW 5 = 0. 

SW 4 = 0. 

SW 5 = 1. 

sat = o. 

BEMFT = EEMFC - EEMFA 
VN 1 = VIF - VEflFC 
VN 2 = VIE - VEMFA 
IM=IMC 
60 CONIINUE 
ENDFBOCEDUEE 

TEBMINAL 

TITLE EASIC DC M01CE SYSTEM 

TIMEE FINIIM = .050. COTBEL = .0000275, FELEL = .00C275,... 

DE1MIN= 1 . 3E- 10 
FINISH IHES'I = 400. C 

PRINT Ml, WlIEFM, I HE ST , 1M, 1MIN I, I MINI, uMlNI ,IAV,1 AV,WMAV, ... 
VAIND, VEIND, VCIND/VN 1 . V N 2 , V AC , V EC , VCC , C 3 , V CD4 , V 0 D 5 , V C E 6 , . 
V5G, VAG, VEG, VCG,VD1D,V02D,VD3C,VC4C,V£5£,VE€D, ... 

IAB, IAEA, ISC, IECA,lCA # ICAA,Ifl,IfcEEG, T ECC N , THCC N 1 

IMA. I MB, IMC, F*E,V1N,VINHAF,IMCAC,VSGF ,1MK,TEM, VKLWK 
LABEL EHASE C CUEEENT - AMES 
OUTPUT TIME, I MC , S k 2 
PAGE XYELCT 
ENE 

BESET EBINT 

LABEL FCWEE SUPPLY CUEEENT TC MCICE - AMPS 
CUTFU1 TIME, 1 M , S « 4 
PAGE XYPLCT 
EEC 

LABEL MCTot SFEZL - K P 1 
OUTPUT TIME, W«6t3,Sik£ 

PAGE XYFLCT 
E NC 

LABEL PCWSk SUPTL* TEEM INAL VCITAGE - VCLTS 
OUTPUT TIME, V3G,S* t 
PAGE XYFLCi 
END 

LABEL TOTAL DEVELOPED TOECUE - C2-IN 
OUTPUT TIME, 1M,SW1 
PAGE XYPLCT 
END 

LA3EL EP-rILTEBEC PHASE C lUSBFAT - AMPS 
OUTPUT TIME, I MC AC , S * 3 
PAGE XYPICI 
EEC 

LABEL VOLTAGE ACEOSS A-C TEBM1KALS - VCLTS 
OUTPUT TIME, V A C , S W 3 
PAGE XYFLCI 
END 

LABEL LP-FILTEEEC FGwDxT SUPPLY CUEEENT -AMPS 
CUTFUT TIME, IM.F,Sk5 
PAGE XYELCI 
EEC 

LAEEL dACK Eilf PHASE C - VCLTS 
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LE; BflTEA 



CSMP 



A 1 



DUTPUT TIKE, VEEFC,SW 1 
PAGE XYFLC1 
END 

LABEL IS CLVELCEEC ICFgUE - CZ-IN 
DUTPUT LIME, TMH,S«3 
PAGE XY£LCT 
END 

-ABEL PC WEE CUTOUT - WATTS 
DUTPUT USE, PHR,SW£ 

PAGE XYEICT 
E NB 

LABEL LP-FILTEEEC SUPPLY TEHMIIsAL VCLIAGE - VCL1S 
OUTPUT TIKE, VSGF,SW1 
PAGE X Y £ LC T 
END 
STCP 

3 JOB 

3ISSFCP. MYDI5S EC * 

DRA W= 1-ENE, 

MODIFY= 1-END ( 3) (SIZE-8.,2.), 

HODIFY = 2-ENC(J) (C V E £ P LOT , CC B. » E E= 0. , 3 . , S I Z E = 8 . , 3- ) , 
SODIF Y = J- END (3) <CiV£F£LOT,CCKNEF = 0.,3.,SIZ£ = e.,3.) 
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Users Guide 



The following procedures are recommended to be used in running the CSMP 
language program given in Appendix B: 

1. Obtain the motor, load and power supply/conditioner data from the 
manufacturer and enter the values in INITIAL section at the beginning 
of the program. 

2. Set parameter KRAMP to the desired input voltage V as follows: 

V 3 

example V = 30 volts, KRAMP = xIO = 15000. Under these conditions, 

a fast terminated-ramp input results. This is equivalent to a step input 
of V volts. 

3. The constant KADJ must be set experimentally to produce the measured motor 
Back EMF voltage. To do this, remove the asterisk in the statement 

* WM = 320.4475 

where 320.4475 is the magnitude of the driver motor speed in rad/sec. (in 
this case, speed = 3060 RPM). Run the program with different values of 
KADJ until the Back EMF across 2 windings agree in magnitude with the 
measured value obtained from the motor. 

4. Set N to the value of the speed reduction from motor shaft to load 
shaft . 

5. Set RS to the value of the power supply interval resistance and LA and 
RA to- the per phase winding values. 

6. Set RSAT and VSAT values for the power transistor "ON" values and RCUT to 
the cut-off resistance value. 



BIO 



7. The value of BM should b« calculated from the no-load data as a 

preliminary value and adjusted experimentally until no-load operation of 
the model agrees with measured no-load values for the motor. Compute BM 
as follows: 

BM = (KT) (KADJ ) (IMNL)/(WMNL) 

For example 

BM = (15.89)(.63)(0.3)/(320.4) 

BM = 0.00937 oz-in/rad per sec. 

Notes: 1. The asterisks in column 1 of the program statement makes the 

statement a comment and disregarded by the CSMP Translator. These statements 
are the result of program development and are included for future application 
in the development of the complete actuator model. 

2. A set of typical output waveforms are included here to indicate 
the plot output capabilities of the program which are in the TERMINAL section 
of the CSMP listing. 

3. The load applied to the motor should be adjusted over the range of 
no-load to near peak load and the performance curves of Figure 12 plotted. 
The CSMP simulation output can then be checked against these results by 
entering BL as follows: 

BL = TL/WM oz-in/rad per sec 

where TL is the load torque (oz-in) WM is the motor speed (rad/sec). 
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POWER OUTOUT - WATTS 




TOTAL DEVELOPED TORQUE - OZ-IN 



□ = PWR 
o = SW5 





MOTOR SPEED - RPM 



□ = TM 
o = SW1 




Figure B1 Step Input Response Waveforms 
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LEGEND 



LEGEND 






LEGEND 






0.6 1 1.6 B 8.6 9 9.6 0 0.85 0.50 0.76 1 1.86 



LP- FILTERED POWER SUPPIY CURRENT -AMPS 




POWER SUPPLY CURRENT TO MOTOR - AMPS 



□ = IMF 
o = SW5 



LEGEND 




PHASE C CURRENT - AMPS 



LEGEND 




Figure B2 Additional Step Response Waveforms 



LEGEND 
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